The signature of dark energy perturbations in galaxy cluster surveys 
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All models of dynamical dark energy possess fluctuations, which affect the number of galaxy 
clusters in the Universe. We have studied the impact of dark energy clustering on the number of 
clusters using a generalization of the spherical collapse model and the Press-Schechter formalism. 
Our statistical analysis is performed in a 7-parameter space using the Fisher matrix method, for 
several hypothetical Sunyaev-Zel'dovich and weak lensing (shear maps) surveys. In some scenarios, 
the impact of these fluctuations is large enough that their effect could already be detected by existing 
instruments such as the South Pole Telescope, when its data is combined with WMAP and SDSS. 
Future observations could go much further and probe the nature of dark energy by distinguishing 
between different models on the basis of their perturbations, not only their expansion histories. 
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Introduction - Over the past few years evidence for the 
accelerated expansion of the Universe has become over- 
whelming: besides supernovas [TJ HJ [31 S], every well- 
measured cosmological observable seems to support the 
conclusion that we no longer live in a decelerated, matter- 
dominated Einstein-de Sitter Universe H [71 HI H] . The 
search for the cause of this acceleration, whether it is yet 
another dark component (dark energy) or if it is due to 
some modification of gravity which makes it repulsive at 
large scales [101 [HI [HI [13] ; has become one of the most 
pressing questions of our times. 

The properties of dark energy, and its effect on cosmo- 
logical observables, are commonly parametrized in terms 
of an equation of state = Pde/Pde with a simple, 
dependence on redshift, Wde{z) = wq -\- Waz/{1 -\- z) 
[14', [15] . Substantial observational efforts are underway 
[16, 111 HI], and many more are being planned [191 120]. 
in order to determine such parameters. 

However, getting the equation of state even with ex- 
cellent precision, in whatever parametrization, would still 
tell us close to nothing about the nature of dark energy 
[21j , or if the acceleration is in fact due to modified grav- 
ity. One way of answering those questions is to look 
for evidence of dark energy perturbations in cosmologi- 
cal observables - evidence which must exist if dark energy 
is anything but the Cosmological Constant [22]. Notice 
that the signal of dark energy perturbations should not 
be confused with the infiuence of smooth dark energy 
on the growth function of dark matter inhomogeneities, 
which is another important, but distinct effect of dark 
energy upon structure formation [331 HH US] • 

Unfortunately, searching for imprints of dark energy 
fluctuations can sometimes be underwhelming: since 
daxk energy dominates only at late times, its impact 
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on the cosmic microwave background (CMB) arises only 
through the integrated Sachs- Wolfe effect, which in turn 
makes the signal both small, and limited to the largest 
scales where cosmic variance defeats the statistical pur- 
pose of the data [26 ] (27 1 (28 1 [29] . 

Dark energy fluctuations also affect the power spec- 
trum P{k) [HI [5S1 [m [301 [31]. In particular, Hu [32] 
used an effective description for the dark energy pres- 
sure perturbations in order to estimate their impact on 
some cosmological observables. In that approximation, 
the pressure perturbation 5pde — ScS^pde, where Soff is 
an effective sound speed squared and 5pde is the energy 
density perturbation of dark energy |69| . This frame- 
work is not ideal, but it is extremely useful if our goal is 
to enlarge the class of models we want to constrain, just 
like any given parametrization of the equation of state 
is not meant to represent actual models, but to simplify 
the task of assessing the power of observations to con- 
strain those models. Along these lines, e.g. Takada ^33] 
regarded Soff as a free parameter and concluded that mea- 
surements of the CMB and power spectrum could detect 
the effects of dark energy inhomogeneities, but only if the 
effective sound speed was very small, which would make 
it behave almost like dust on the scales probed by the 
CMB and P{k). 

That leaves the nonlinear regime of structure forma- 
tion (collapsed structures, or halos) as one of the few 
venues left where we could expect the influence of dark 
energy to be measurable. As the initially overdense 
regions attract the surrounding matter and experience 
gravitational collapse, the steep gravitational potentials 
at the center of the dark matter halos deform the distri- 
bution of dark energy. And since pressure is a key ingre- 
dient in halo formation, dark energy could have an en- 
hanced effect on the masses, times of formation and rela- 
tive abundances of halos ,34, 35, 36, 37, 38 , 39 , 40 , 41, 4^ 
as well as voids [43l [4^. Notice that both P{k) and 
galaxy cluster counts are potentially sensitive tests of 
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dark energy perturbations. However, while the former 
tests the long-range, infra-red properties of the underly- 
ing theory of dark energy in the linear regime, the latter 
is a test of the small-scale, ultraviolet, nonlinear regime 
of the same unknown theory. They are, therefore, com- 
plementary rather than competing tests of the nature of 
dark energy - see also [35^ on using cluster counts to 
constrain models of modified gravity. 

Methodology ~ Counts of galaxy clusters have long been 
regarded as a crucial cosmological observable ^\ ST] , in 
particular for purposes of constraining dark energy mod- 
els L4_8, 49, 50, ^ Ea ESI EH- In the present work we 
study the sensitivity of cluster counts to dark energy per- 
turbations. We employ a set of seven cosmological pa- 
rameters that we allow to vary: h (the value of the Hub- 
ble expansion rate, where Hq = 100 h km Mpc~^), 
flm (the amount of matter relative to the critical den- 
sity), flh (density ratio of baryons), as (rms matter fluc- 
tuations on scales of 8 h"'^ Mpc), wq, Wa and Scff (we 
assume flatness throughout.) The abundance of halos is 
then estimated with the Press- Schechter formalism [55] . 
using the top-hat spherical collapse model 155] to com- 
pute the critical density at the time of collapse - see also 
[SI |46l EH EH] , and, in particular, |40l [42] . 

Our statistical analysis relies on the Fisher matrix, 
which was computed in this 7-dimensional parameter 
space for several hypothetical Sunyaev-Zel'dovich (SZ) 
and weak lensing (WL) surveys of galaxy clusters, and 
a suitable flducial dynamical dark energy (DDE) model 
with wq = —1.1 and Wa = 0.5. We have checked that the 
impact of dark energy clustering scales like 1 1 -I- | , which 
is natural since in the ACDM limit there are no dark en- 
ergy perturbations. The other flducial values we used 
are h = 0.72, fit = 0.05, erg = 0.76 and fl^ = 0.25. This 
model lies near the best-flt region in the parameter space 
of several recent joint analyses of cosmological observa- 
tions [59l[60]. Notice that, even though our DDE model 
happens to experience "phantom crossing" (w — —1) at 
a redshift z = 0.25, there is no associated instability in 
the spherical collapse model, hence the sensitivity to dark 
energy is not artiflcially enhanced in this scenario com- 
pared to similar scenarios without phantom crossing - 
see, for instance, [HT]. Finally, for completeness we stud- 
ied three scenarios of DDE, with flducial values Soff — 0, 
SoS = 0.5 and Soff = —0.75. 

The main limitation to the sensitivity of number counts 
is the halo mass below which clusters cannot be detected 
with the given instrument and strategy. For the i-th bin 
the number of observed clusters is given in terms of the 
mass function as: 

where AH. is the solid angle subtended by the survey. The 
limiting mass Mmin is approximately constant in redshift 
for surveys which employ the SZ effect, but in WL cluster 
surveys the limiting mass grows with redshift due to the 




Figure 1: Joint Constraints on s^ff and f2m from near-future 
(nf) surveys. The top panels show marginalized 68% C.L. 
limits from SZ-type surveys, and the bottom panels show the 
limits from WL-type surveys, for the fiducial models with 
Se// = (left panels) and Seff ~ —0.75 (right panels.) The 
thick solid (black in color version) lines correspond to cluster 
counts limits (68% C.L., marginalized); the thick dashed (red) 
lines assume an additional weak prior on Seff, f^s^ff = 1; the 
dark (blue) ellipses correspond to adding the "COSMO" set of 
priors (see text) to the cluster Fisher matrix; and the smaller, 
light (green) ellipses correspond to adding both the weak prior 
on Seff and the COSMO priors. 



declining number of background galaxies |621 163j . 

We consider three hypothetical scenarios for these two 
types of surveys: near-future (nf) ones, such as the 
SPT/DES dS], which are supposed to cover about 4000 
deg^; future (f) ones, such as the LSST [18], which will 
cover about half the sky (18.000 deg^); and far-future (ff) 
ones, such as EUCLID/JDEM [IHIIMI, which will cover 
up to 75% of the sky. 

We have deflned limiting masses by a simple two- 
parameter function: log;^g Mijiin/(/i^^M0) = X + Pz. We 
have set X^f = 14.2, = 14.1 and Xff = 13.7, and 
= in all cases - which means we assume exactly 
flat selection functions for the SZ-like surveys, clearly a 
very rough approximation. For the WL surveys we have 
set X^f^ = 14.0, Xf^ = 13.5 and Ajf^ = 13.2, and 
pWL Q Q jj-^ ^-|^-|^ cases. With these limiting masses, 
we obtained (for ACDM) total numbers of clusters of: 
7000 and 4600 respectively for the near-future SZ and 
WL surveys; ~ 280.000 for future surveys (both SZ and 
WL); and ^ 1.5 x 10^ for far-future surveys (both SZ and 
WL). 

We have used different binnings in mass and in red- 
shift in each case: for near-future SZ and WL surveys 
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Figure 2: Joint Constraints on Se// and from WL surveys 
in the fiducial model with s^ff — 0, for future surveys (upper 
panel) and far-future surveys (lower.) The colors are the same 
as in Fig. 1. 



we used 3 mass bins, and 10 redshift bins up to z — 1.5; 
for future SZ and WL surveys we used 5 mass bins and 
15 redshift bins up to z = 1.5; and for the far-future 
SZ and WL surveys we used 8 mass bins and 25 red- 
shift bins up to 2; = 2.0. Our results do not change 
significantly if fewer mass bins are used, and we chose 
fat redshift bins in order to minimize sample ( "2-halo" ) 
variance [63] |65l [66]. In fact, we only consider Poisson 
(shot) noise in the covariance of the number of detected 
halos: {{Ni - N,){Nj - Nj)) = SijNi, where Ni (Ni) is 
the actual (expected) number of halos in the i-th bin. 



The Fisher matrix is then given by: 

i 

where 9°' are any one of the 7 cosmological parameters. If 
all other parameters are fixed, the Icr (68% C.L.) limit on 
the parameter 9'^ is given by \/l/Faa- If the other param- 
eters are marginalized over (integrated out of the total 
PDF), then the Icr hmit on 61" is given by yJ{F-^)aa ■ 

We have also considered the effect of adding an exter- 
nal ("COSMO") set of cosmological priors: WMAP 5y 
(through the shift parameter R) [67 , BAO [6 , HST [BH] 
and a baryon fraction n^h? = 0.02273 ± 0.00062, also 
from WMAP 5y data [67]. 

Results and Discussion - If the fiducial model is ACDM, 
then naturally the sensitivity to perturbations (and to 
the sound speed) of dark energy vanishes. For the dy- 
namical dark energy fiducial model we find that the im- 
pact of the dark energy perturbations is dramatically dif- 
ferent depending on the value of the dark energy pressure 
perturbation. If the pressure is positive, zero or negative, 
then the impact on the mass function (and on the num- 
ber counts) is respectively very small, small or very large 
(see Figs. 1-3.) This means that the likelihood function 
is more peaked for smaller values of the sound speed. 

For positive dark energy pressure, our results are sim- 
ilar to what would be expected from "quintessence" 
(canonical scalar field) models of dark energy, since in 
that case the Jeans length is approximately equal to 
the Compton wavelength of the dark energy field, Aj ~ 
Cs/wde (cs — y^Seff is the sound speed.) In that case, 
dark energy clusters only on horizon scales, and its effect 
on cluster counts is minute. 

For near-vanishing dark energy pressure (sound speed 
very small) our results are similar to power-spectrum con- 
straints like those by obtained by, e.g., Takada [33]. How- 
ever, because the typical Jeans lengths are still typically 
much larger than cluster scales, the limits from number 
counts are weaker than the limits from P{k). 

For negative pressures (se// < 0) the sensitivity of 
galaxy cluster counts to dark energy perturbations in- 
creases dramatically. In fact, this would be true for the 
sensitivity of any large-scale structure observable (e.g., 
CMB and power spectrum), if it were not for the fact 
that it is probably nonsensical to consider negative sound 
speeds squared at the linear level. However, in collapsed 
regions and small scales there is nothing that forbids 
negative pressures ^34l HI] ~ in fact, the bag model of 
hadronic physics is one microphysical example that sup- 
ports the use and interpretation of negative pressure in 
cosmology. 

Our main conclusions are: 

1. Ongoing surveys such as the South Pole 
Telescope, which will detect a few thousands 
of clusters, are already sensitive enough to 
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Figure 3: Joint Constraints on wo and Wa from near-future 
SZ surveys, for the fiducial models with Se// = (left) and 
Seff = —0.75 (right.) Legends for the thick solid and dashed 
lines, as well as for the ellipses, are identical to the previous 
figures. In addition, computing the 68% C.L. limits without 
the dark energy clustering parameter s^ff leads to either the 
thin dashed (black) or dot-dashed (brown) lines, whether or 
not one includes the COSMO set of priors. 




Figure 4: Joint Constraints on Q,n and ag from SZ (top) and 
WL (bottom), near-future (left) and future (right) surveys. 
In all panels the fiducial model has Se// = 0. Legends are 
identical to Fig. 3 

detect dark energy perturbations (in the ab- 
sence of strong systematics) if Seff < 0. 
Dark energy perturbations also weaken the 



detection of the other dark energy parame- 
ters determined by the SPT, such as wq and 
Wa'- in particular, the marginalized figure-of- 
merit {1/ \/AwqAw^) is substantially smaller 
50%) compared to the case where dark en- 
ergy perturbations are neglected. This points 
to a degeneracy between the effective sound 
speed and the equation of state which is not 
unexpected, since one determines the pres- 
sure on small scales, while the other sets the 
pressure on large scales. 

2. Future surveys, such as the LSST, may 
or may not find evidence of dark energy per- 
turbations. But even in the worst scenario 
{seff > 0, weakest dependence on the sound 
speed) the constraints on the main cluster pa- 
rameters (such as flrn and erg) are substan- 
tially weakened by this additional "nuisance 
parameter" . Furthermore, the marginalized 
limits on the equation of state parameters wq 
and Wa are affected at the level of 10-30%, 
depending on the scenario - see Figs. 3-4. 

3. Surveys in the more distant future, such 
as ESA's Euclid or NASA's JDEM, will de- 
tect the imprints of the fluctuations of dark 
energy, but in the most conservative scenario 
the constraints would only be strong enough 
to determine whether the dark energy pres- 
sure in collapsed regions is positive or nega- 
tive. 

Therefore, if future data points to anything other than 
a Cosmological Constant, cluster surveys will be a key 
observation to determine the nature of dark energy. But 
even if cluster counts turn out to impose weak constraints 
on dark energy clustering, this would still be a minor 
ingredient affecting the halo mass function. Our results 
show that the addition variance (nuisance) introduced by 
dark energy clustering in the determination of the other 
cosmological parameters through galaxy cluster counts 
cannot be dismissed. 
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